Quantum dots with two electrons: Singlet-triplet transitions 
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The magnetic character of the ground-state of two electrons on a double quantum dot, connected 
in series to left and right single-channel leads, is considered. By solving exactly for the spectrum of 
the two interacting electrons, it is found that the coupling to the continuum of propagating states 
on the leads, in conjunction with the electron-electron interactions, may result in a delocalization of 
the bound state of the two electrons. This, in turn, reduces significantly the range of the Coulomb 
interaction parameters over which singlet-triplet transitions can be realized. It is also found that 
the coupling to the leads favors the singlet ground-state. 



I. INTRODUCTION 



The spin state of the many electron ground-state is 
determined by the interplay between the kinetic and the 
electrostatic interactions. According to Hund's law the 
ground-state of electrons in a partially filled shell of an 
atom has the maximal possible spin, in order to mini- 
mize the electrostatic repulsion. On the other hand, An- 
derson's super-exchange antiferromagnetic interaction, 
which favors zero total spin, arises from the reduction 
in the ground-state energy brought about by the hop- 
ping of electrons between adjacent ions. Another exam- 
ple is realized in gases made of homonuclear diatomic 
molecules, in which the distance between the nuclei deter- 
mines the magnetic properties of the gas; in other words, 
the ground-state of the molecule can be either a triplet 
or a singlet. 

The possibility to study these rules as function of con- 
trolled parameters, and in particular to observe devia- 
tions from them, has been opened up in recent years with 
the intensive experimental and theoretical investigations 
of small confined systems. Devices based on quantum 
dots formed in GaAs heterostructures allow to probe the 
electronic states in situations where there are only a few 
electrons in the system, as function of, e.g., the number of 
electrons in the sample, the spacing of the single-electron 
energy levels, the internal structure of the dots (e. g. 
the distance and the coupling between internal parts) 
or the coupling of the system to external leads. Thus, 
modulating the single-particle spectrum by a magnetic 
fielcL.|applied perpendicularly to the plane of the electron 
gas,ElH has produced a structure of the conductance peaks 
that has been interpreted in terms of singlet-Jjciplet tran- 
sitions of the last two electrons in the dotSu (see also 
Rcf. |@). This has been observed on vertical quantum 
dotsp and also in lateral ones.cl The ground-state spin 
of chaotic quantum dots has been studied by tracing the 



conductance peak spacing as function of a weak parallel 
magnetic field (which couples primarily to the spins). Q 
A numerical investigation of such a configuration ana- 
lyzed the influence of the exchange interaction on the 
peak structurc.El In a similar fashion, the singict-triplet 
transitions in such dots have been attriubted to avoided- 
crossings in the many-electron states, and the relation 
between those -and kinks in the conductance pattern has 
been explored B 

When the device consists of two or more coupled dots, 
another controllable parameter conies into the play: the 
interdot coupling. Experimentally,li3 it has been found 
that this coupling shows up in the conductance peak posi- 
tions. Theoretically, the effect of the interdot distance on 
the states of a few electrons confined in a parabolic paten- 
tial has been analyzead-, using mean-field methods,EJ tha 
Kohn-Sham equationJ13 and numerical diagonalization.ll3 
Different spin states of double-quantum-dots systems 
have been also studied using the numerical rcnormaliza- 
tion group method £j 

Here we present an exact analytical solution for the 
eigenvalues and eigenfunctions of two electrons on a 
double-dot system, which is coupled to two single-channel 
leads. The two electrons interact while they are on 
the dots, and we include in the calculation the direct 
Coulomb and the exchange interactions between the dots, 
and the on-site Hubbard interactions on the dots. We 
obtain the condition for the singlet-triplet transitions of 
the ground-state energies, and in particular examine the 
role of the delocalization effect of jLheJnteractions: It has 
been fouadjn our previous works JIj^Ej as well as in other 
studies,llZrlla that the interplay between the hybridization 
of the localized single-particle states on the dot with the 
propagating states on the leads, and the electron-electron 
interactions, may promote one of the electrons to the con- 
tinuum. If the system had only one doubly bound state, 
then the above 'promotion' results in a 'delocalization' of 
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the ground state. Once delocalized, the singlet and the 
triplet states become degenerate. 

Our earlier work included only the case with on-site 
Coulomb repulsion. In that case, the bound ground state 
is a singlet. As the coupling of the dots with the leads was 
increased, the parameter range where this bound state 
exists was decreased, and the delocalization was accom- 
panied by a transition from the singlet to a degenerate 
singlet-triplet 'metallic' state. In the present paper we 
add the effects of direct and exchange interactions. In the. 
absence of the coupling to the leads, it is well knownta 
that the ground state changes from a singlet to a triplet 
upon increasing the exchange interaction. Here we inves- 
tigate what happens to that singlet-triplet transition in 
the presence of the coupling to the leads. 

After discussing our model Hamiltonian in Sec. II, 
we present the exact two-electron solution in Sec. III. 
Section IV then discusses the singlet-triplet transition, 
mainly for a special choice of the parameters where it is 
easiest to explore the solution analytically. Our results 
are then summarized in Sec. V. 



II. HAMILTONIAN 

The Hamiltonian of the model reads 

H=Hsp+K. (1) 

The single-particle tight binding Hamiltonian is 

T-isp = J2 '^i4aCia - ^ijcLcja, (2) 
ia ijo- 

where c}^ creates an electron with spin a on site i. The 
site energies ei are different from zero only on the 'quan- 
tum dots', which can be viewed as 'impurities' on the lat- 
tice. In what follows, we specifically consider two neigh- 
boring quantum dots, which will be denoted by £ and 
r, and study the symmetric case where ei^r — co- The 
hopping matrix elements are divided into three kinds: 
the hopping among the 'dots' is denoted by tjj (= t^), 
the hopping between a dot and a neighboring site on the 
'lead' is denoted by to (e. g. tn for i 7^ r), and the hop- 
ping between neighboring 'lead' sites is denoted by i = 1, 
setting the units of energy. 

The Coulomb interactions are assumed to exist only 
among electrons which sit on the dots. Generally, this 
interaction has the form 

= Tji C;g.c|j^/CniCT'Cno-, (3) 

ijmn (7(7' 

where 

ri-" = JdvJ dr'T;(r-r')^r(r)^*(r')^m(r')¥'n(r), (4) 

while ifi{r) is the Wannier wave function of an elec- 
tron localized at site i. It is usually assumed that the 



dominant terms-will be those in which the indices are 
equal pair-wise.L3 Here we follow j-the parametrization 
used in magnetic coupling studies ,Eil in which the case 
i = j = m = n is treated separately. Neglecting the coef- 
ficient with i j 7^ m = n, which (when negative) leads 
to the superconductivity vertex, we are left with three 
possible parameters: 

2r|| - u, 

2r|! = v, 2r;] = K, i^j, (5) 

where U, V and K are the (intra-dot) Hubbard, (inter- 
dot) direct, and (inter-dot) exchange interactions, respec- 
tively. Then 

Hint = Uy^ni|?T,i| 
i 

-|-i(2V-K)^ninj-K^Si-Sj, (6) 

i#j i#j 

where 

a 

Si = y^Cj^gg-g-'Cicr' (7) 
0-0-' 

and Gaa' is the vector of Pauli's spin matrices. It is im- 
portant to note that with the above approximation for 
the Coulomb interaction vertex, the interaction Hamil- 
tonian becomes spin-dependent. 

A similar Hamiltonian has been used to study the ef- 
fects of the spin states on the condjicLiag properties of 
confined mesoscopic quantum dots.aQ'Q'cl In these stud- 
ies, i represents some orbital state on the dot. In Ref. 

it has been argued that when the Thouless conduc- 
tance of the confined (mesoscopic) system is large, the 
Coulomb vertices can be parametrized in terms of two 
coupling constants, independent of the orbital indices. In 
terms of our parameters, these are the charging energy, 
Ec = (2V — K)/4, and the exchange energy, J = K, while 
U = V -|- K. It is straightforward to apply our results to 
that case. 

In our special case of the double dot, the energies F in 
Eq. (^ differ from zero only for i, ] = (, r, and thus 

Hint = \J{h^nii + flr^hrl) 

+ (V - ^)nihr - 2KS, • Sr. (8) 
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III. THE TWO-ELECTRON 'MOLECULE' 

We will confine ourselves to a double-dot system (a 
'molecule'), containing two electrons. The system is mod- 
eled by two identical single- level impurities, each having 
an on-site single-particle energy level ei = er — eo- The 
two impurities are coupled to one another by the inter- 
dot matrix element to- When the molecule is isolated 
{to — 0), there are three degenerate triplet states, of en- 
ergy 2eo -I- V — K, and three singlet states. One of the 
latter has the energy 2eo-|-U, while the other two energies 
aral3 



/U- V-K\2 



(9) 



Examination of the interaction matrix elements, Eqs. (Q) 
and (|), shows that U>V>OandU>K. Mor cover, 
the exchange interaction K involves the square of the 
overlap matrix element of the two impurities, and hence 
is of the order of t^. Therefore one concludes that the 
lowest singlet energy is given by the minus sign in (|^). 
The triplet states energy is below the lowest singlet one 
provided that 



2tl < K(U- V + K). 



(10) 



One notes that when the direct and exchange Coulomb 
vertices V and K are disregarded, the energy of the sin- 
glet state is lowered by the interdot kinetic energy, pro- 
ducing the, Anderson super-exchange antiferromagnetic 
couplingjHa ^ 4t|)/U. The triplet state becomes the 
ground-state cpce the exchange energy K wins over this 
kinetic energycSIa Finally, with the choice of parameters 
in which U = V -f K (see above) , the condition becomes 
tf) < (see also Ref. ||). Both to and K decay expo- 
nentially with the interdot distance. However, Eq. (Q) 
indicates that K is roughly of order t^, and thus decays 
faster. Therefore, one might expect a transition from 
the triplet to the singlet ground state as this distance 
increases. 

The main purpose of the present paper is to study how 
the criterion for the singlet-triplet transitions, Eq. (10), 
is modified when the double quantum dot is connected to 
external leads, on which it is assumed that there are no 
electron-electron interactions. In other words, we study 
the changes in these transitions due to coupling the dot 
to a continuum of propagating states. A similar prob- 
lem was discussed in Ref. ^ in the framework of the 
mean-field approximation, within a single-lead geometry. 
Here we first derive exactly the ground-state energy of two 
electrons, and then compare the singlet and the triplet 
ground-state energies. 

For further simplification, we describe the external 
leads by single channel one dimensional chains, with 
nearest,, neighbor hopping t = I. In our previous 
work,Ej'E3 we have shown that the spectrum and the wave 
functions of two interacting electrons can be obtained in 



terms of the energy spectrum and the wave functions of 
the single-particle Hamiltonian. For our model, cf. Eq. 
(^), the latter can be divided into even and odd solu- 
tions, denoted by "e" and "o". Of particular interest 
here will be the regions in the parameter plane {eo — 7}, 
where 7 = ig, where the spectrum has bound states, see 
Fig. ^ For eo < 7 — 2-|-t_D, there appears one (even) 
bound (single-particle) state below the band of propa- 
gating states, of energy eg < —2; For even smaller values 
of eo, such that cq < j — 2 — tu there appears a second 
(odd) bound state below the band, of energy e^, with 



2-7 



I-7V 2 / 1- 
^0 = ^0 — to, = eo 



■7 

to 



7-1, 



(11) 



Similarly, for eo > 2 — j — to there appears the first 
bound state above the band, while for eo > 2 — -f + tjj 
there are two bound states above the band. Clearly, a 
necessary condition to have the two electrons bound in a 
triplet state is the existence of two distinct single-electron 
bound states, so that each electron occupies a different 
'state'. For simplicity, we shall confine ourselves to the 
case in which both occur below the continuum, i.e., to 
the lowest region II in Fig. |l|. 




FIG. 1. The single-particle bound energies. The roman 
numbers denote the number of bound states. Here ti3=0.4. 

It has been shown in Refs. |l^ and |l^ that when inter- 
actions between the two electrons are allowed, they may 
lead to the delocalization of one (or both) of the electrons 
from the doubly bound states. Thus, for example ,Ej when 
V=K=0, in part of region I in Fig. one finds that the 
doubly bound ground state is replaced by a ground state 
in which one of the electrons is shifted to the band. In 
these parameter regions there are no two-electron bound 
states, and the system may be called 'metallic'. Indeed, 
such interaction-induced delocalization efi^ects have been 
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recently observed in quantum dots containing two po- 
tential minimaEj This effect has also been found for two 
interacting electrons moving in a one-dimensional peri- 
odic structure, EZI |and for spinless fermions on strongly 
disordered chains Ji3 In contrast to region I, in region II 
there potentially exist two doubly bound states, i. e. the 
singlet and the triplet, and our goal is to study which of 
these states is the ground state in the presence of this 
interaction-induced delocalization. 

To solve for the spectrum of the two interacting elec- 
trons we proceed as follows. Let us denote the eigenstates 
and the eigenvalues of the single-particle Hamiltonian ^ 
by (pain) and ea, respectively, where n is the site index. 
Consider first the two-electron states with spin Sz — ±1, 



ab 

ab nn' 

^\[\Y.^nn'<^nJn'M^ '^^il, (12) 

nn' 

where the amplitudes X^„, — — X^/„ are antisymmetric 
in the coordinates (the origin is half way between the 
dots), and X^rm' l^mi'P = 1- The Schrodinger equation 
for the two electrons then yields 



= (V - K) 



ah 



rAnWdn'). (13) 



Hence, the two-particle energies, E, are given by the so- 
lutions of the equation 



1 



V-K 2 



ab 



E - ea - efc 



(14) 



Making now use of the symmetry properties of 0a, i-e., 
C(0 = and 0°(r) = -(/)°(^), we see that the 

only contributions to the sum come from the cases in 
which a is even, and b is odd, and vice versa. For the 
sake of brevity, in the following expressions we replace 

0a(^)l = l</^a(^)l (* = 6,0) by I'/'al- Consequently, we 
may write the result in the form 



1 



V-K 



Goo(E), 



(15) 



where Geo is a noninteracting two-electron Green's func- 
tion. 



Gco(E) = 4 



2EE 

i,j=e,o k 



E - e^, - efc 



E 



kk- 



f E - efc - efc/ 



(16) 



(In writing down this equation, we have assumed the 
existence of two single-particle bound states, as men- 
tioned above.) Here the subscript B denotes a single- 
particle bound state, and k refers to a band state, with 
Ck — — 2cosfc. In the continuum of the band energies, 
there is no need to distinguish between e^ and e^. Also, 
the sum over all states k is divided into the sum over the 
even propagating states, and the sum over the odd ones. 
We give in the Appendix the explicit expressions for the 
eigenstates required for the calculation of G. 

Equation ( p^ ) is an implicit equation for the exact two- 
electrons eigenenergies E. We will postpone the discus- 
sion of these solutions, and consider now the two-electron 
states with Sz = 0, 



ab 

b nn' 

x^^/c|^^c]^, I io), 



ab 



(17) 



with J2nn' l^nn'P — 1- Eor the triplet state, the ampli- 
tudes X^, are antisymmetric in the site indices. Then 
(inserting Eq. ( pT[ ) in the Schrodinger equation) the en- 
ergies are again given by Eqs. (|l5|) and (|l^). For the 
singlet states, X°„, — X^,„ and the Schrodinger equa- 
tion yields 



x^ = E uE^' 



(l>a{i)(l>b{i) 



ab \ i—i,r 

E - ea - eb 



E - ea - eb 

KinWbin') 



(18) 



We use this equation for n, n' = £, £, n, n! = r, r, and 



n, ri! — i?, r, and find two families of singlet solutions: (i) 
XO = -XO^, XO = 0, for which 



— Gco(E); 



(ii) XO^ = XO^, with 

1 - (U + V + K)(Gec(E) + Goo(E)) 
+ 4U(V + K)Gec(E)Goo(E) =0. 



(19) 



(20) 



[Note that the last equation includes twice the number 
of solutions as Eqs. ( [l5| ) and (p^.j Here, Goc,oo are non- 
interacting, two-particle Green's functions, which consist 
of the even and odd (with respect to interchanging the 
dots) solutions of the single-particle spectrum, respec- 
tively, with 



G'oc(oo)(E) — 2 



E-2e 



= (o) 
B 



fc 



E 



o(o) 



e/c 



Ae(o)|2| i,o(o)|2 



2f- E-efe-efc, 
kk' 



(21) 
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We next determine the ground-state energy, starting 
with the triplet states, whose energies are given by Eq. 
(p^). We assume that we are in the lower region II of 
Fig. 1^, where there are two bound states below the 
band. (Note that in the regions marked and I in 
Fig. 1^, where there is at most one single-particle bound 
state, the triplet 'bound' state will always lie in the con- 
tinuum.) The function Goo(E), Eq. ([l^), has the fol- 
lowing behavior As E approaches — oo, it goes to 
zero from below. At E= -|- Cg, it diverges to — oo, 
jumps to -foo as E crosses that value, and then de- 
creases, as E approaches the bottom of the two-electron 
continuum states,|4pcated at — 2 -I- e|. As discussed in 
our earlier workllslla, in the thermodynamic limit of infi- 
nite 'leads' Geo has a finite value at this band thresh- 
old, due to the ^-dependence of (jfi}"^ at the impuri- 
ties. The triplet ground-state energy is where Geo crosses 
(V — K)~^. Hence, there will be a two-electron bound 
state only when (V - K)"^ < Geo(-2 -I- Cq). It follows 
that there are values of the direct and exchange Coulomb 
couplings such that the ground triplet state is not boimd, 
but lies in the continuum. 

We now turn to the singlet states, again assuming the 
existence of both and eg . Consider first the solutions 
given by Eq. ([l9| ) . Since U> V (and V-K) , the lowest so- 
lution of this equation lies above the lowest solution of the 
triplet state, (which is given by the same function Geo)- 
Hence, we need not consider anymore the states given by 
(|l9|). To explore the other family of singlet solutions, it 
is convenient to re-writc Eq. (EOh in the form 



Similarly to the behavior of Geo(E), the right-hand-side 
of this equation starts at very small negative values when 
E tends to — oo. It then diverges to — oo as E approaches 
2eg, jumps to -|-oo as E crosses that value, diverges again 
to — oo as E— > 263, then jumps to -l-oo, and finally de- 
creases towards a finite value as E approaches the bot- 
tom of the two-electron continuum. It follows that Eq. 
( p3| ) has always a bound energy solution. Moreover, if 
Goo + Goo is negative at E= Et = -I- e^, then that 
solution Eg lies below Ex, i.e., the ground-state is a sin- 
glet. This is indeed the case, as is shown in the Appendix 
[Eq. ([aT^)]. This is in accordance with the general rule, 
which states that in order for the ground-state to be a 
triplet, the exchange Coulomb energy has to overcome 
the kinetic energy. 

For the sake of clarity of the presentation, we will carry 
the rest of the analysis to lowest order in the coupling 
to the leads, 7. In that case, it is possible to derive sim- 
ple ex press ions for the two-particle Green's functions, see 
Eqs. (A13). Using those equations, we find that the sin- 
glet energies are given by 



E-2eo+7(e-"'' +e-"°) 



U + V + K 



U- V-K\2 



) +4iD7(e-"° -e-"°), (24) 



and the triplet energies are given by 

E - 2eo + 7(e""' + e""") = V - K, 



(25) 



1 / 1 



1 



U + V + K 



4 \ Goo(E) Goo(E) 



1 / 1 



\ 2 / \ 2 

1 \ /U-V-K\ 



\ 16 \^Goo(E) Goo(E)^ 



(22) 



Since the behavior of Goo and Goo as function of E is 
similar to that of Goo described above, it follows from 
Eq. (^2|) that the lowest singlet state energy obeys that 
equation with the minus sign. 



IV. SINGLET-TRIPLET TRANSITIONS 

In order to decide when the lowest bound state of the 
two electrons is a singlet or a triplet, we need to (i) deter- 
mine for which values of the Coulomb parameters Eqs. 
( p^ ) and (HI) have bound solutions; and (ii) to compare 
these two solutions, when they exist. Consider as an ex- 
ample the case in which there is only the on-site Hubbard 
interaction, that is, V=K=0. Then the triplet bound 
state has the energy Et = Cg + ^Iq- (p^. The 

singlet energy. Eg, in that case is given by the lowest 
solution of 



^ = Goo(E) 



Goo(E). 



(23) 



where a'''" is related to the corresponding E (with indices 
S or T) via 



E - e;:^° = -2cosha°'°. 



(26) 



Let us examine the case in which in the absence of the 
coupling to the leads, the singlet and the triplet ground- 
state energies are equal, i.e., 4t|) = 2K(U-V+K), see Eq. 
(0). Then 



Eg - Et = 7(e""T + e" 
4tD 



U - V + 3K 



"T _ e-"s _ e-"s) 
7(e-"°s -e-"s). (27) 



The last term on the right-hand-side of this equation is 
negative. This follows from Eq. ( p6| ) and the fact that 
£3 < £%■ fo'" the first term, we use again Eq. (|2^), to 
write it in the form (Et - Eg)7((e2"s - l)-i + (e^^s - 
1)"^). Hence, Eg < Et and the singlet is prefered. 

The above discussion shows that the coupling to the 
continuum of propagating states enhances the tendency 
of the two electrons to form a singlet state, in the situ- 
ation where in the absence of that coupling, the singlet 
and the triplet states are degenerate. In order to inves- 
tigate whether this tendency persists for other choices 
of parameters (and at the same time to keep the calcu- 
lations tractable) we will now confine ourselves to the 
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choice U=V+K. In this case, again to leading order in 7, 
the equation for the singlet energies (24) reads 



/s(E) = V + K, 

/s(E) = E-2eo + 2^z5 + 27e- 



(28) 



Similarly, Eq. ( |25|) for the triplet energies can be written 
as 



/t(E)=V-K, 

/t(E) = E-2eo+7(e-«°+e-"°). 



(29) 



Let us first determine for which parameters these equa- 
tions yield bound, two-electron energies. To this end, we 
consider Eqs. ( p^ ) and ( p9| ) at the bottom of the two- 
electron continuum, E=— 2 + eg. The first of these equa- 
tions will have a bound state for /s(— 2 + eg) > V + K; 
the second will have such a solution when /t(— 2 + eg) > 
V — K. These conditions are plotted in Fig. |^ as the 
thick lines there. A bound triplet state exists in regions 
I-I-III, below the heavy line of positive slope. A singlet 
bound state exists in I-I-II, below the heavy line of nega- 
tive slope. In region IV, there are no bound states; both 
the triplet and the singlet states are in the continuum, 
and their energy is about the same. Crossing the line 
between regions II and IV (or III and IV) thus corre- 
sponds to the delocalization transition discussed above, 
from a singlet (triplet) bound ground state to a degener- 
ate 'metallic' state. This transition is the most striking 
effect of the coupling to the leads, and we expect it to ap- 
pear irrespective of the quantitative approximation used 
in Fig. |. 

In region I one has to compare the singlet energy with 
the triplet one. These two become degenerate along the 
diamond curve, whose equation is derived from (26), (28) 
and (|2~ 



2eo-(l-7)AB, 



V-K = -^(B-1) + ^(B + 1) 

A= -(K-ic), 

7 



B = Jl + 4/[A(A + e° -6| 



(30) 



This line is almost vertical, with Es < Et to its left, and 
Es > Ft to its right. The conclusion is that, as long as 
there exists a bound state of the two electrons, then this 
line moves slightly to the right as 7 is increased from zero 
(when the line was at K— tjj). However, the coupling to 
the continuum states delocalizes the electrons, making 
the two states degenerate over a large part of the param- 
eter plane {V — K}, i. e. region IV in Fig. ^. With all 
other parameters fixed, one might expect that increasing 
the distance between the dots causes a decrease in to, 
in V and in K, thus causing a shift towards to lower left 
side of Fig. ||, towards a bound singlet ground state. 




FIG. 2. The phase diagram of the two-electron bound 
states, for to = 0.4, eo = —2.6, and 7 = 0.1, see text. 



V. DISCUSSION 



We have derived analytical expressions for the spec- 
trum of two interacting electrons on a simplified model 
for a double quantum-dot. When the dot is decoupled 
from the external leads, it is straightforward to obtain 
this spectrum, and discuss the criterion for its ground- 
state to be a singlet or a triplet. The question we have ad- 
dressed is how this criterion is modified when the single- 
particle states which are localized on the dot are coupled 
to the continuum of extended states on the leads. A typ- 
ical example of our results is shown in Fig. |^: as long 
as the the electron-electron interactions do not delocalize 
the ground-state, then the location for the singlet-triplet 
transition shifts continuously with the coupling to the 
leads, 7. In that case, one can still say that the singlet 
state is the ground one provided that the kinetic energy 
dominates over the exchange energy. In a way, the cou- 
pling to the leads enhances the kinetic energy, and there- 
fore it slightly favors the singlet ground state. However, 
the coupling to the leads has a much more drastic effect: 
over a significant part of the parameter space (which con- 
sists of the Coulomb couplings, the single-particle ener- 
gies on the dots, and the coupling to the leads), e. g. 
region IV in Fig. |^, the interplay between the coupling 
to the leads and the electron-electron interactions delo- 
calizes one or both electrons. Then, the singlet and the 
triplet states are degenerate. In such a situation, the 
bound state (which may be either a singlet or a triplet) 
disappears and the ground state becomes a degenerate 
singlet-triplet 'metallic' state. We believe that this delo- 
calization effect should be taken into consideration in the 
analyzes of experimental data related to this question. 
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APPENDIX A: THE TWO-ELECTRON GREEN'S 
FUNCTIONS 

We first list the eigenfunctions required for the calcu- 
lation of the nonintcracting, two-particle Green's func- 
tions. These are needed only on the double quantum 
dot, that is, on the sites r and £. Writing the even and 
odd bound-state energies, Eqs. (O), in the form 



~2coshK'' 



the wave functions on the dot sites is 



2(1+ e2«'= '' - 1 



(Al) 



(A2) 



The band states have been calculated assuming periodic 
boundary conditions for a system of N sites (these include 
the leads). Then 



,0,0 1 2 



7 sin k 



(A3) 



with efe = — 2cosA;, for both the even and odd states. In 
calculating the sums over k in the two-particle Green's 
functions, we shall use the continuum limit, dividing the 
^-integrations on the even (odd) functions by 2. 

We next derive Goo, Goo, and Goo- It is convenient 
first to calculate the function 



Writing 



and using Eq. 



oj = — 2cosha, 
we find 



LO < -2. 



(A4) 



(A5) 



(e"+'=°'° - 1 + 7) (6"+'^"'° - 1) ' 
where we have also used 



(A6) 



(A7) 



Exploiting this result, wc now introduce the function _F, 



F°-°(c.) = |0-°|2 + (c. 



_ 1 1 
^ 2 ^ 1 + (7/(6"+'*'' 



(A8) 



Note that 



1 



(7/(e 



(A9) 



It is now straightforward to show, using Eqs. (pi]), (p^), 
and (^ that 



Goo(E) - 2 



(f -(E - 6° ))- 
E-2e?, 



E 



fefc 



- (E - efe - efcO(E - e% - ek){E ^ e% ~ ek') ' 



(AlO) 



with an analogous result for Goo, with e replaced by o, 
and 



Geo(E) 



E 



f -(E - e^)F°(E - e|) 
rk\^m^{e°B-ek){e%-ek,) 



kk 



efc)(E 



ek') 



(All) 



For energies E below the bottom of the two-electron con- 
tinuum, i.e., E< —2 -|- eg, the double sum on k and k' in 



(AlC) is negative. Using Eq. (A9) for the first terms in 
the equations for Geo and Goo, it follows that 



Gco(eB + eB) + Goo(eB + eB) <0. 



(A12) 



This result is used to show that when the only Coulomb 
coupling is the Hubbard U, the singlet is always the 
ground state. 

Up to this point, the results were given for general 7. 
To lowest order in the cou pling to th e lead s, 7, we may 
discard the double sums in ( A10| ) and ( All ). Then, using 
|), we find 



1 



Geo(E) 
1 

2Goc(E)' 
1 



2Goo(E) 
where we have used Eq. (& 



E-2eo + 7(e~" + e"" 
E-2eo + 2ti5 + 27e~"°, 



(A13) 



* Permanent address: School of Physics and Astronomy, 
Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel. 
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